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In an open system, the geometric phase should be described by a distribution. We show that a geometric 
phase distribution for open system dynamics is in general ambiguous, but the imposition of reasonable physical 
constraints on the environment and its coupling with the system yields a unique geometric phase distribution 
that applies even for mixed states, non-unitary dynamics, and non-cyclic evolutions. 
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Introduction. — The geometric phase (GP) identi- 
fies the portion of an overall (abelian or non-abelian yj) 
phase shift of a quantum state that is due to the path of the 
state through projective Hilbert space; the GP and the dy- 
namic phase combine to give the aggregate phase shift of the 
state that may be inferred by interferometric or other phase- 
sensitive methods. GP theory has been rigorously formulated 
for the general case of non-adiabatic |4], non-cyclic JpL and 
non-unitary evolution (without quantum jumps) |0, 0, 0] of 
a pure state, but the importance of GP in realistic systems, 
for example in the context of adiabatic quantum computation 
has motivated recent important research into GP in 
open systems 1 11, RTlHOlJi', 19, 20]. Quantum jump 
(or trajectory) analyses have been applied to certain physical 
systems, which show how the GP for a closed system can be 
modified under open system dynamics iiH Eill . and a rigor- 
ous Kraus operator approach has been applied to define GP 
for general open system evolution 1 23] . These studies note the 
importance of GP beyond closed-system, unitary evolution of 
pure states. 

In this paper we argue that a complete description of abelian 
GP in open systems has to identify the appropriate measure of 
phase distribution. We develop a theory of GP distributions 
for mixed states, non-unitary dynamics, and non-cyclic evo- 
lutions. We show that, without further contraints, the GP dis- 
tribution is ambiguous: an operational definition of GP that 
would resolve this ambiguity is not attainable because the GP 
is a non-linear functional of the state. The ambiguity of the 
phase distribution is rather subtle: we show that previous def- 
initions of GP distributions and its spread implicitly assume a 
particular form of phase distribution. The imposition of rea- 
sonable physical contraints on the environment and its cou- 
pling with the system yields a unique GP distribution by tak- 
ing the decomposition of the density matrix |23] into account. 

Definition of geometric phase distributions for an open 
system. — Interferometric or other phase-sensitive measure- 
ments allow inference of the phase shift of a state, but sepa- 
rating geometric and dynamic components of the phase is not 
straightforward. For a pure state \'4>{t)) — V{t)\ilj{Q)) E Ti 
which is propagated by an arbitrary time-dependent evolution 
operator V{t) (not necessarily cyclic or unitary), the mathe- 
matical definition of the geometric phase functional /3[i/'] is 



given by 0] 

which is meaningful only for Z[ijj] ^ 0. The dynamic phase 
functional of 

removes the dynamic phase from the total phase shift asso- 
ciated with V{t). In some cases, the dynamic phase can be 
eliminated via interferometry of a state that follows a super- 
position of two paths, with the dynamic phase along each path 
being the additive inverse of the other |24, 25], but this can- 
cellation of dynamic phase is not always achievable. 

Non-unitary evolution of the type V{t) may not satisfy the 
axioms of completely positive (CP) maps, which guarantee 
that a positive-definite operator on Hilbert space such as the 
density operator p is mapped to a positive definite operator 
with identical trace, and linearity is preserved. Thus the GP 
should be established for general CP maps, not just for non- 
unitary evolution llsUldllTll . A physical picture for the CP 
emerges by considering a system S with Hilbert space Tis 
and a reservoir R (or set of ancillae) with Hilbert space TYr, 
and joint Hilbert space Ti = Tis ® TCr. At some initial time 
t — the joint density factorizes, psr(O) = ps(0) (8) Pr(0), 
and the unitary operator C/sr(^) of the system-nreservoir (Sh-R) 
is generated by a hamiltonian i/sR(0 such that ps,R{t) = 
UsR{t)psR{0)Uln{t). Dynamics for S alone is obtained by 
tracing over R; i.e. 

Ps(i)-TrR(c/sRpsR(0){/sV) (3) 

is a CP mapping, which can be decomposed into a sum of 
mappings corresponding to various measurement records ob- 
tained by readouts of the environments. For a particular initial 
state ps(0) = Y.s9s\^s){'4's\ andpR(O) = this 
Kraus decomposition corresponds to 

ps(i) = Pr'i^^r,bAt)m{^s\Ul^At) (4) 

with Kraus operators Ub^^r = (&r|C^sr(OI'')' where {|6r)} is 
some basis for the Hilbert space of the reservoir. This cor- 
responds to an incoherent mixture of non-unitarily evolving 
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states \ipbR,r.s) = C4g 1 V's), weighted with the initial prob- 
abilities PrQs- For each individual state the GP is given by 
Eq. Q with Z[^|Jb^^r,s] oc {tpbR,r,s{0)\fpbR,r,s{t)). Because 
lV'fcR,r,s(0)) = {bR\r)\tps), it is obvious that Z[?/;6R,r,s] = 
whenever the basis state |6r) is orthogonal to the reservoir 
state \r), so that the associated GP is not well defined. To 
avoid this problem we start by choosing a different basis set 
{|^r(^))} for each term in the sum over r in Eq. (0}, such 
that {bfi{r)\r) = 6^^^. Consequently, when one introduces a 
distribution for complex numbers of the form 

P(z) = J2 ^('^' " ^[^bR,r,s]) , 2 G C , (5) 

bR,r,s 

with w{r, s, 6r) being arbitrary weight functions, all terms 
with 5r 7^ do not contribute to any moment (z") = 
/ z''P{z)d^z and can therefore be omitted. Hence, it is suf- 
ficient to keep only terms with 6r = so that the complex 
number distribution induced by Eq. becomes 

Pz{z) = Y,PrQsSiz-Z[lbr,s]), (6) 

r,s 

with \ijjr,sit)) = {r\UsR{t)\r) lips). For a total system Sh-R 
with a continuous spectrum the sums in this distribution would 
be replaced by integrals. 

The definition for a corresponding GP distribution faces the 
same subtleties that arise for any phase distribution. Usually, 
a positive operator-valued measure (POVM) is introduced to 
describe a phase distribution for quantum systems. However, 
through the dynamic phase functional D [^] , the GP depends 
non-linearly on the states of the system, so that it is generally 
not possible to provide a GP POVM based on linear opera- 
tors. Instead, one has to construct a GP phase distribution 
differently. We investigate two natural definitions of GP dis- 
tributions. The first possibility to introduce a GP distribution 
is to derive the GP directly from the complex number distri- 
bution Pz of Eq. (|6ji; for example the mean GP is given by the 
first moment i26ll 

l(^>z| \J2r,sPr'lsZ[^pr,s]\ ' 

By definition, £)[V'r.s] becomes unity if each \ipr.s) is paral- 
lel transported. In this case definition Q can be rewritten 
as exp(i(/3)) = Tr(C/sR(t)psR(0))/|Tr([/sR(t)psR(0))|, which 
coincides with the definition given in Ref. 113]. 

A second possibility to introduce a GP distribution is mo- 
tivated by Holevo's approach to moments of a phase distribu- 
tion fr?] . While the definition (Q depends on the modulus of 
Z[i(j], too, this is not necessarily desirable. Instead, one can 
introduce phase distribution of the form 

The corresponding first moment is given by 

(e''^)^(e-)ff =^p,g,e''3['^-=l. (9) 

r,s 



An advantage of this expression is that its phase can be consid- 
ered as a mean GP, while its modulus is related to a measure 

W ^ \(e''^)\-^ - 1 (10) 

for the spread of the GP. Eqs. Q and (|5Jl yield different results 
for the mean GP, which reflects the choice available in obtain- 
ing an average for the GP from the same distribution. Ideally, 
in each run one of the states \tpr,s) is realized and leads to a 
well defined complex value for Z[tpr,s]- It is then only a mat- 
ter of definition how the average of the complex values over 
all runs is performed. 

Initially pure systems and density matrix decomposition. — 
To concentrate on the effects of the reservoir we focus on an 
initially pure system with ps(0) — |'0s)(V's|i so that 

Pz{z)=Y.Pr5{z-Z[i^r]), (11) 
r 

with \il]r{t)) = {r\U5R{t)\r)\tps) . An immediate consequence 
for a reservoir in a pure state is that the GP distribution is sharp 
(i.e. a 5 function). An important example of the reservoir be- 
ing initially in a pure state is any energy state, in particular its 
ground state. 

For the reservoir in a mixed state, however, the non-linear 
dependence of the dynamic phase functional D[tp] on |?/;) 
leads to an ambiguity in Eqs. ([^ and (|8} if the density matrix 
can be decomposed into mixtures of two different sets of states 
1 23]. The GP then not only depends on the choice of distribu- 
tion but also on the density matrix decomposition. We argue 
here that, with respect to the mean GP, the introduction of 
physical constraints can resolve both ambiguities. Naturally- 
occurring reservoirs do not exhibit coherence between differ- 
ent energy levels, so we assume the density is block diagonal 
in the energy basis (with block sizes determined by degenera- 
cies); the thermal reservoir with density operator cxp{—/3Hr) 
is a typical example of such a reservoir density matrix. It is 
then physically reasonable to only admit decompositions of 
the density matrix which differ with respect to the decompo- 
sition in degenerate subspaces. As the dynamic phase func- 
tional D{E) is identical for all states sharing the same energy 
eigenvalue E, the contribution of the respective subspace He 
takes the form 

Pr5{z - i?(i;)(V's|(r|C/sR|r)|Vs)) ■ (12) 

reUm 

Hence, the corresponding contribution to (z)z can be written 
as Trs.-Hg (C/srp(0)), which is independent of the decompo- 
sition. This is not the case for the mean GP in the Holevo 
measure, so that the resolution of the decomposition problem 
favors the choice of the measure Pz(^) for a GP distribution. 
We note that for higher moments the decomposition ambigu- 
ity remains for both choices of GP distribution. However, we 
will see below that for a weakly coupled reservoir the GP dis- 
tributions coincide and are independent of the decomposition. 
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Explicit expression for a weakly coupled reservoir — To 
derive an explicit expression for GP when the system's evo- 
lution is described by a family of CP maps, we start with a 
total hamiltonian of the form H{t) = Hs{t) + Hj!_ + Hi with 
constant weak coupling Hi and time-dependent system hamil- 
tonian Hs{t), which generates a unitary evolution UsR{t) = 
rexp(— i dt' H{t')) of the total system. The reservoir is 
initially in a mixture of eigenstates |r) of the time-independent 
hamiltonian Hr. In the interaction picture, Usn{t) = 
Us{t)UK{t)U{t), with Unit) = exp(-itiJR) and Us{t) ^ 
rexp(— i dt' Hs{t')). Standard second-order perturbation 
theory then leads to U{t) = 1 + A + B + 0{Hf) with A = 
-i/o dt'iJi(t') as well as B ee -/^ dt'/J' dt" Hi{t')Hi{t") 
and the interaction-picture hamiltonian Hi = uIu^HiUsUr. 

Instead of explicitly deriving closed expressions for the op- 
erators A and B, we relate them to a corresponding CP map- 
ping of the reduced density matrix ps, which is assumed to be 
of Lindbladtype fH, 

Ps = -i[Hs+AH, ps]-^iLlL^ps+psLiL^-2L^psLi) 

a 

(13) 

where are the jump operators. The hermitean operator 
AH describes any energy shifts (such as the Lamb shift for 
atom-light interaction) associated with Sh-R interaction. On 
the other hand, we can also directly calculate ps to second 
order in Hi by use of Eq. 0. An elementary calculation leads 
to 

PS = -i[Hs,Ps] + {Us{B)i,ulps (14) 
+UsY,Pr{My)UipsUs{r'\A^\r)ul+H.c.} . 

r,r' 

Expectation values are denoted by {■ ■ ■)q for a state l^Jq); 
in particular, (• • • )s refers to (■0s | • ■ ■ IV's) and (• • • )r to 
TrR(/9R(0) • • • ). Comparing Eq. (I14> with the Lindblad form 
(I13> allows us to identify 

Us{B)rUI = -iAH - LiL^ (15) 

a 

while the operator A is related to the jump terms in a non- 
trivial way. 

To keep the presentation concise we now focus on a cou- 
pling of the form Hi = — Rfj.Sfi, where i?^ and S*^ are 
operators which act only on Hr and TCs, respectively. Fur- 
thermore, we assume that 

(r|i?^|r) =OVr , (16) 

which is the case for energy-transfering Sh-R interactions, for 
instance. An immediate consequence is that (A),. = 0. It 
follows that Z[ijjr] = Z[xPs]{l + (AZ)^,s) with 

AZ^-^-^—^+i fdt'{B^AHs+AHsB) , (17) 



with Hs. = U^HsUs and AHs = Hs-{Hs)s- Using Eq. O 
one easily finds the following expression for the moments as- 
sociated with each of the two different GP distributions intro- 
duced above, 

= = e'"'^['^^l(l +inIm(AZ>,,,(o)) . (18) 

Result ( I18> has some interesting consequences. First, Pz and 
Pn generate exactly the same moments so that (e''^) — e'^'^^ to 
second order in Hi. Hence for a weakly coupled reservoir Pz 
and Pii are equivalent and, since I (e''3) I « 1, the GP distribu- 
tion is sharp. Second, as a direct consequence of Eq. ilbi . Pz 
and Pii do not depend on the operator A, so that the GP distri- 
bution does not depend on the jump operators appearing in the 
master equation. Third, as a consequence of the "linearity" of 
expressions to lowest non-trivial order in perturbation theory, 
the moments dlSt are invariant under a change of the reser- 
voir basis. Hence, they are independent of the decomposition 
ofpR(O). 

Explicit Calculations of Berry phase. — As illustration we 
consider the GP distribution defined above for some specific 
physical examples. First, we discuss the case of a two-level 
atom with ground state \g) and excited state |e) interact- 
ing with a thermal radiation reservoir. The corresponding 
hamiltonian is given by Hs = — (huj / 2) {\e) {e\ — \g){g\), 
and the jump operators of Eq. (I13> are given by Li ~ 
\/jo(ji + T}\g) {e\ and L2 = ^7on|e) {g\. Here, 70 denotes 
the spontaneous emission rate and n the thermal mean num- 
ber of resonant photons 1 29j] . At temperature T = we have 
n = so that this model reduces to a description of sponta- 
neous emission in vacuum. The operator B introduced above 
then reads B = — 7o(|e)(e| + nl). For simplicity we have 
omitted the Lamb shift f3(J\. For an initial state of the form 
IV's) = cos I |e) + sin | \g) we find for the (sharp) GP at time 
t = 2tt/u! the expression 

{(3) = (3[^s]+7r^-smH , (19) 

UJ 

where the GP for a closed system is given by /3['0s] = 
27rsin^(0/2). An interesting feature of this result is that the 
mean GP does not depend on the temperature, even though 
for T > the radiation reservoir is in a mixed state. This is 
a consequence of the n-dependence of B being proportional 
to the identity: the effect of thermal fluctuations, which in- 
duce incoherent absorption and emission of thermal photons 
at equal rates, do exactly cancel each other The asymmetric 
effect of spontaneous emission, however, leads to a change in 
the GP 

Alternatively, one can calculate the GP by solving Eq. (I13> 
in terms of Kraus operators Ki, without explicit reference to 
the reservoir, which provides an exact calculation of the GP 
distribution. This simply amounts to seeking for operators Ki 
and probability weights pi for which J2iPi^iPs{^)K} is a 
solution of the master equation |31]. For the present case, the 
Kraus operators are given by 
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ifo = e-¥*|.g)(.g|+e't*-7.*|e)(e| (20) 

= v/l-e-27„t|g)(e| (21) 

= e-**-''"*|.g)(5H-e'^*|e)(e| (22) 

i^a = Vl-e-27„t|e)(g| (23) 

with7„ = (2n+ 1)70 and weights po = Pi = (n+l)/(2n+l) 
as well as p2 — P3 — n/ {2n +1). The operators and 
are related to the jump operators and, because of Ki{Q) = 
K3 (0) = 0, do not contribute to the GP. Substituting Kq and 
K2 into Eqs. and (|9} we calculate the GP distribution at 
time t = 27r/w to be 

P„W^P..(e"-i^)+«.(e'--^) (24, 



s (e 



For zero temperature we have 71 = and therefore p2 = 0. 
Both expressions then predict a sharp GP 



(/3>='r+— ln(^s|e- 
270 



-27r7ocr2/w 



(25) 



This result agrees with the expression found in Ref. |21] and, 
to first order in 70, also with the result il9\ based on the 
weakly coupled reservoir Any difference between expi(/3) 
and (cxp i/3) is of second order in 70. Also for finite tempera- 
tures, the two exact results still agree with the weak coupling 
result jl9t to first order in 70. Hence, any dependence on the 
temperature through n is of higher order in 70. 

Another illustrative case is that of phase damping which 
can be described by a jump operator of the form Li = 
-ya(|e)(e| — \g){g\) (and therefore B cx 1), where a denotes 
the phase damping rate. This jump operator can be derived 
from a coupling to a non-resonant reservoir of harmonic os- 
cillators with effective interaction hamiltonian Hi ~ SqRo ~ 
Gz gi'A'^i' where is the annihilation operator of the ith 
oscillator and gi the corresponding effective coupling param- 
eter In thermal equilibrium we have (R) — g.i{a\ai) ^ 
so that condition M6\ is violated. Consequently the (trivial) 
result predicted by Eq. ( I18> is spurious. 

We can again compare this to an exact calculation based on 
the Kraus operators 
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'\9){9\ 



|g)(5|+re'^*|e)(e| (26) 

ie-^ 
r 



\e){e\ 



(27) 



with r = (1 + \/l — exp(— 2Q;t))^/2 ^nd weights pQ — pi — 
1/2. For brevity we only discuss the first moments which are 
given by 



,i(/3> 



3uj 



■ cos t/ sm 



sin^6'(i cos sin^ 



(28) 



(29) 



for t = 2tt/uj. In contrast to Eq. (II 8> . these moments include 
non-trivial corrections and differ from each other as well as 
from the result of Ref. |l2l|. By Eq. (O they indicate a GP 
spread of « IGtt^ sin'' 6a/ (9lj) for phase damping. 

In summary we have established a theory for GP distribu- 
tions based on operational considerations that reduces to the 
result of Samuel and Bhandari |5] for no-jump non-unitary 
evolution and employs a Kraus operator analysis. We ad- 
dress and resolve ambiguities concerning decomposition of 
the density matrix and GP by incorporating reasonable as- 
sumptions about the reservoir and solve specifically for spon- 
taneous emission and phase damping of a two-level atom. The 
resultant GP distributions are derived and the ambiguities over 
phase mean and spread are consequences of choices of how to 
define a phase distribution. The underlying method for defin- 
ing the GP distribution is now clear, but of course choices 
of the distributions are ultimately determined by experimental 
considerations. Future work will consider the GP distribution 
for non-abelian GP. 
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